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We examine a synthetic way of constructing the grain size distribution in the interstellar medium (ISM). First 
we formulate a synthetic grain size distribution composed of three grain size distributions processed with the 
following mechanisms that govern the grain size distribution in the Milky Way: (i) grain growth by accretion 
and coagulation in dense clouds, (ii) supernova shock destruction by sputtering in diffuse ISM, and (iii) shattering 
driven by turbulence in diffuse ISM. Then, we examine if the observational grain size distribution in the Milky 
Way (called MRN) is successfully synthesized or not. We find that the three components actually synthesize the 
MRN grain size distribution in the sense that the deficiency of small grains by (i) and (ii) is compensated by the 
production of small grains by (iii). The fraction of each contribution to the total grain processing of (i), (ii), and (iii) 
(i.e., the relative importance of the three contributions to all grain processing mechanisms) is 30-50%, 2Q-AQ%, 
and 10^-0%, respectively. We also show that the Milky Way extinction curve is reproduced with the synthetic grain 
size distributions. 
Key words: cosmic dust, interstellar medium, grain size distribution, extinction, Milky Way. 



1. Introduction 

Dust grains are important in some physical processes in 
the interstellar medium (ISM). For example, they dominate 
the absorption and scattering of the stellar light, affecting the 
radiative transfer in the ISM. The extinction (absorption + 
scattering) by dust in the ISM as a function of wavelength is 
called extinction curve (Wickramasinghe, 1967; Hoyle and 
Wickramasinghe, 1991; Draine, 2003 forreview). Extinction 
curves are important not only in basic radiative processes in 
the ISM but also in interpreting observational data: part of 
stellar light in a galaxy is scattered or absorbed by dust grains 
within the galaxy in a wavelength-dependent way according 
to the extinction curve. Therefore, to derive the intrinsic stel- 
lar spectral energy distribution of a galaxy, we always have 
to correct for dust extinction by considering the extinction 
curve (Calzetti, 2001). 

Extinction curves generally reflect the grain composition 
and the grain size distribution. Mathis et al. (1977, hereafter 
MRN) show that a mixture of silicate and graphite dust, as 
originally proposed by Hoyle and Wickramasinghe (1969), 
with a grain size distribution (number of grains per grain 
radius) proportional to a -35 , where a is the grain radius 
(a ~ 0.001-0.25 /im), reproduces the Milky Way extinc- 
tion curve. Pei (1992) shows that the extinction curves in the 
Magellanic Clouds are also explained by the same power-law 
grain size distribution (i.e., oc a~ 35 ) with different abun- 
dance ratios between silicate and graphite. Kim et al. (1994) 
and Weingartner and Draine (2001) have applied more de- 
tailed fit to the Milky Way extinction curve in order to obtain 
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the grain size distribution. Although their grain size distri- 
butions deviate from the MRN size distribution, the overall 
trend from small to large grain sizes roughly follows a power 
law with an index near to —3.5. Therefore, the MRN grain 
size distribution is still valid as a first approximation of the 
interstellar grain size distribution in the Milky Way. 

What regulates or determines the grain size distribution? 
There are some possible processes that actively and rapidly 
modify the grain size distribution. Hellyer (1970) shows that 
the collisional fragmentation of dust grains finally leads to 
a power-law grain size distribution similar to the MRN size 
distribution (see also Bishop and Searle, 1983). In fact, Hi- 
rashita and Yan (2009) show that such a fragmentation and 
disruption process (or shattering) can be driven efficiently 
by turbulence in the diffuse ISM. However, they also show 
that grain velocities are strongly dependent on grain size; 
as a result, the grain size distribution does not converge to 
a simple power-law after shattering. Moreover, dust grains 
are also processed by other mechanisms. Various authors 
show that the increase of dust mass in the Milky Way ISM 
is mainly governed by grain growth through the accretion 
of gas phase metals onto the grains (we call elements com- 
posing dust grains "metals"). (Dwek, 1998; Inoue, 2003; 
Zhukovskaef al., 2008; Draine, 2009; Inoue, 2011; Asano et 
al, 2012). In the dense ISM, coagulation also occurs, mak- 
ing the grain sizes larger (e.g., Hirashita and Yan, 2009). In 
the diffuse ISM phase, interstellar shocks associated with su- 
pernova (SN) remnants (simply called SN shocks in this pa- 
per) destroy dust grains, especially small ones, by sputtering 
(e.g., McKee, 1989). Shattering also occurs in SN shocks 
(Jones etal., 1996). 

Modeling the evolution of the grain size distribution in the 
ISM is a challenging problem because a variety of processes 
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are concerned as mentioned above. Those processes are also 
related to the multi-phase nature of the ISM. Liffman and 
Clayton (1989) calculate the evolution of grain size distri- 
butions by taking into account grain growth and shock de- 
struction. However, their method could not treat disruptive 
and coagulative processes (i.e., shattering and coagulation). 
O'Donnell and Mathis (1997) also model the evolution of 
grain size distribution in a multi-phase ISM, taking into ac- 
count shattering and coagulation in addition to the processes 
considered in Liffman and Clayton (1989). They use the ex- 
tinction curve and the depletion of gas-phase metals as quan- 
tities to be compared with observations. Although their mod- 
els are broadly successful, the fit to the ultraviolet extinction 
curve is poor, which they attribute to the errors caused by 
their adopted optical constants. They also show that inclu- 
sion of molecular clouds in addition to diffuse ISM phases 
improves the fit to the observed depletion, but they did not 
explicitly show the effects of molecular clouds on the extinc- 
tion curve. Yamasawa et al. (201 1) have recently calculated 
the evolution of the grain size distribution in the early stage 
of galaxy evolution by considering the ejection of dust from 
SNe and subsequent destruction in SN shocks. Since they 
focus on the early stage, they did not include other processes 
such as grain growth and disruption (shattering), which are 
important in solar-metallicity environments such as in the 
Milky Way (Hirashita and Yan, 2009). 

Comparing theoretical grain size distributions with obser- 
vations is not a trivial procedure. In a line of sight, we always 
observe a mixture of grains processed in various ISM phases. 
Therefore, a "synthetic" grain size distribution, which is 
made by summing typical grain size distributions in indi- 
vidual ISM phases with certain weights, is to be compared 
with observations. In this paper, we first formulate a syn- 
thetic way of reproducing the grain size distribution. Then, 
we carry out a fitting of synthetic grain size distributions to 
the observational grain size distribution, in order to obtain 
the relative importance of individual grain processing mech- 
anisms. We do not model the multi-phase ISM in detail, but 
our fitting contains the information on the weights (i.e., rel- 
ative importance) of different grain processing mechanisms, 
which depend on the ISM phase. 

This paper is organized as follows. In Section [2] we 
explain our synthetic method to reconstruct the grain size 
distribution in the ISM. In Section [3] we fit our synthetic 
models to the observational grain size distribution in the 
Milky Way and examine if the fitting is successful or not. 
In Section |U after we discuss our results, we calculate the 
extinction curves to examine if our synthesized grain size 
distributions are consistent with the observed extinction or 
not. In Section[5] we give our conclusions. 

2. Synthetic grain size distribution 

As explained in Introduction, we "synthesize" the obser- 
vational grain size distribution (here, the MRN size distri- 
bution) by summing some representative grain size distribu- 
tions in various ISM phases. These representative grain size 
distributions are explained in Section |2~TI 

First, the ISM is divided into two parts: one is the 
part where the grain processing is occurring (called "grain - 
processing region"), and the other is the area where the 



grains already processed in the various grain-processing re- 
gions are well mixed (called "mixing region"). The mass 
fractions of the former and the latter regions are, respec- 
tively, /proc and 1 — iproc- It is reasonable to assume that 
the grain size distribution in the mixing region should be the 
mean grain size distribution in the ISM (Section l2~2l . 

We assume that all grains are spherical with material 
density s; thus, the grain mass m is expressed as m — 
^ira 3 s. Although coagulation may produce porous grains 
(e.g., Ormel et al, 2009), we neglect the effects of porosity 
and assume all grains to be compact. Two grain species are 
treated in this paper; silicate and graphite. To avoid com- 
plexity arising from compound species, we treat these two 
species separately. This separate treatment is also practi- 
cal in this paper as we (and other authors usually) assume 
that the observed extinction curve can be fitted with the two 
species (Section l4~4l ). Before being processed, the grain size 
distribution is assumed to be MRN: a power-law function 
with power index —r(r = 3.5), and upper and lower bounds 
for the grain radii (whose values are determined below) a m i n 
and a max , respectively: 
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nMRN(fl) = 0. The grain size distribution is defined so that 
7t,mrn(o) da is the number of grains whose sizes are between 
a and a + da per hydrogen nucleus. The dust mass density, 
Pd, is related to the metallicity Z (the mass fraction of el- 
ements heavier than helium in the ISM) and the hydrogen 
number density njj as (Hirashita and Kuo, 2011) 



»-H?<h){l 



n-R, 



(2) 



where mx is the atomic mass of the key element X (X = Si 
for silicate and C for graphite), fx is the mass fraction of 
X in the dust, £ is the fraction of element X in gas phase 
(i.e., the fraction 1 — £ is in dust phase), and (X/H) is the 
solar abundance relative to hydrogen in number density. The 
metallicity is assumed to be solar (Z = Z Q ). 

We fix the maximum grain radius as a max = 0.25 /jm 
(MRN). Although the lower bound of the grain size is 
poorly determined from the extinction curve (Weingartner 
and Draine, 2001), we assume that a m i n = 0.3 nm, since 
a large number of very small grains are indeed necessary 
to explain the mid-infrared excess of the dust emission in 
the Milky Way (Draine and Li, 2001). For the other pa- 
rameters, we follow Hirashita (2012). We assume that 0.75 
of Si is condensed into silicate (i.e., £ = 0.25) while 0.85 
(i.e., £ = 0.15) of C is included into graphite. Those 
values are roughly consistent with the observed depletion 
(e.g., Savage and Sembach, 1996), and reproduce the Milky 
Way extinction curve (Section 14.4b . We adopt the follow- 
ing abundances for Si and C: (Si/H)© = 3.55 x 10 -5 and 
(C/H) Q = 3.63 x 10~ 4 . We assume that Si occupies a mass 
fraction of 0.166 (fx = 0.166) in silicate while C is the only 
element composing graphite (fx = 1). We adopt s = 3.3 
and 2.26 g cm~ 3 for silicate and graphite, respectively. 

By using the MRN size distribution as the initial condi- 
tion, we calculate the evolution of grain size distribution by 
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the various processes treated in Section 12.11 In the numeri- 
cal calculation, the grains going out of the radius range be- 
tween a m i n and a max are removed from the calculation (the 
removed mass fraction is < 1%). The processes considered 
are (i) "grain growth" - grain growth by accretion and coag- 
ulation in dense medium, (ii) "shock destruction" - destruc- 
tion by sputtering in SN shocks, and (iii) "grain disruption" - 
grain disruption by shattering in interstellar turbulence . Shat- 
tering in SN shocks (Jones et al., 1996) could be included as 
a separate component, but in our framework, it is not pos- 
sible to separately constrain the contributions from the two 
shattering mechanisms because both shattering mechanisms 
(turbulence and SN shocks) selectively destroy grains with 
a > 0.03 /im and increase smaller grains, predicting sim- 
ilar grain size distributions. Thus, we simply assume that 
the size distribution of shattered grains, whatever the shatter- 
ing mechanism may be, is represented by the one adopted in 
Section l2~T31 

Although our fitting procedures are based on grain size 
distributions, we should keep in mind that observational con- 
straints on the grain size distribution is mainly obtained by 
extinction curves. Weingartner & Draine (2001) performed 
a detailed fit to the Milky Way extinction curve. However, 
the grain size distributions derived by them broadly follow 
an MRN-like power law, although there are bumps and dips 
at some sizes. We will examine the consistency with the ex- 
tinction curve later in Section H~4l 
2.1 Processes considered 

2.1.1 Grain growth Grain growth occurs in the dense 
ISM, especially in molecular clouds, through the accretion 
of metals (called accretion) and the sticking of grains (called 
coagulation). The change of grain size distribution by grain 
growth has been considered in our previous paper (Hirashita, 
2012). The grain size distribution after grain growth is de- 
noted as n gIOW (a, t gvow ), where i grow is the duration of grain 
growth. We adopt t gIOW = 10 and 30 Myr based on typ- 
ical lifetime of molecular clouds (e.g., Lada et al., 2010). 
The metallicity in the Milky Way is high enough to allow 
complete depletion of grain-composing materials onto dust 
grains in ~ 10 Myr. Thus, the total masses of silicate and 
graphite become 1.33 (= 1/0.75) and 1.18 (= 1/0.85) times 
as large as the initial values, respectively (recall that the dust 
mass becomes 1/(1 — £) times as much if all the dust grains 
accrete all the gas-phase metals). The difference in the grain 
size distribution between t gIOW = 10 and 30 Myr is predom- 
inantly caused by coagulation rather than accretion. 

2.1.2 Shock destruction We calculate the change of 
grain size distribution by SN shock destruction in a medium 
swept by a SN shock, following Nozawa et al. (2006). All 
SN explosions are represented by an explosion of a star 
which has a mass of 20 M Q at the zero-age main sequence, 
and the SN explosion energy is assumed to be 10 51 erg. For 
the ISM, we adopt a hydrogen number density of 0.3 cm -3 
(since the destruction is predominant in the diffuse ISM; Mc- 
Kee, 1989), and solar metallicity. The calculation of grain 
destruction is performed until the shock velocity is deceler- 
ated down to 100 km s _1 (8 x 10 4 yr after the explosion). We 
apply the material properties of Mg2Si04 and carbonaceous 
dust in Nozawa et al. (2006) for silicate and graphite, re- 
spectively. We denote the grain size distribution after shock 



destruction by n s hock(a)- The destroyed mass fractions of 
silicate and graphite are 0.38 and 0.27, respectively. 

2.1.3 Disruption Grain motions driven by interstellar 
turbulence lead to grain disruption (shattering) in the diffuse 
ISM (Yan et al., 2004; Hirashita and Yan, 2009). Among 
the various ISM phases, dust grains can acquire the largest 
velocity dispersion in a warm ionized medium (WIM). We 
recalculated the results of earlier workers based on our as- 
sumed initial conditions. We adopt the same grain velocity 
dispersions and hydrogen number density (nn = 0.1 cm~ 3 ) 
in the WIM as adopted in Hirashita and Yan (2009). The 
fragments are assumed to follow a power-law size distribu- 
tion with a power index of —3.3 (Jones et al., 1996; Hirashita 
and Yan, 2009). We denote the grain size distribution after 
disruption as ndisr(a, idisr), where idisr is the duration of 
shattering in the WIM. The lifetime of WIM is estimated to 
be a few Myr from the recombination timescale and the life- 
time of ionizing stars (Hirashita and Yan, 2009). Thus, we 
adopt tdisr = 3 and 10 Myr for our calculation in causing 
moderate and significant disruption. 
2.2 Synthesizing the grain size distribution 

In the beginning of this section, we introduced the mass 
fraction (/ pl0 c) of ISM hosting grains which are now being 
processed ("grain-processing region"). The mean grain size 
distribution over all the grain-processing region, n synt (a), 
can be synthesized with the processed grain size distribu- 
tions, n grow (a, ig r ow) (grain size distribution after grain 
growth with a growth duration of t gmw ), n s hock(a) (grain 
size distribution after shock destruction), and ndisr(a, ^disr) 
(grain size distribution after disruption with a shattering du- 
ration Of tdisr): 

"syntv^v — / grow^grow (**, £grow) ~r /disr^disr (0t) ^disr) 

+/shock^shock(a), (3) 

where f gIO w, /disr and /shock are the mass fractions of 
medium hosting, respectively, grain growth, disruption, and 
grain destruction in the grain-processing region. We call 
^synt(a) "synthetic grain size distribution". If both species 
are spatially well mixed, they would have common values 

f° r /growj /shocks an d Jdisr- 

The mean grain size distribution in the ISM is denoted as 
n mea „(a) and expressed as 

"mcan(^/ — v / proc)^'incan(^) T Jproc^synt {^/j \ v 

since it is assumed that the grain size distribution in the mix- 
ing region has already become the mean grain size distri- 
bution. By assumption, the mean size distribution is MRN: 
n mean (a) ~ nMRN(a)- This condition is equivalent to 



n syn t{a) = "mrn (a). 



(5) 



In the Milky Way ISM, since the grain mass is roughly in 
equilibrium between the growth in clouds and the destruction 
by SN shocks (Inoue, 201 1), we apply f glow Ri = /shock Ri, 
where R\ is the fraction of dust mass growth in clouds 
(0.33 and 0.18 for silicate and graphite, respectively; Sec- 
tion [2HTTJ, and i?2 is the destroyed fraction of dust in a SN 
blast (0.38 and 0.27 for silicate and graphite, respectively; 
Section l2.1.2b . Thus, we put a constraint, 



/shock = (Rl/R2)fg 



(6) 
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We approximately adopt R1/R2 = 0.8 as a mean value 
between silicate and graphite. As mentioned above, if the 
two species (silicate and graphite) are spatially well mixed, 
both species would have common values for f grow , /shock, 
and /disr- Thus, we adopt a single value for i?i/i?2- 
Using the above constraints, Eq. (0 is reduced to 



^synt 



(a) = /growls (a) + /disr«-disr(a, *disr), (7) 



where %, s (a) = n grow (a, t grow ) + (i?i/i? 2 )^ s hock(a). 
Thus, we treat / gr0 w and /disr as free parameters. We de- 
fine the sum of all the fractions as 



/tot — /grow i /shock T" /di: 

R±' 
R 2 , 



1 



/grow ~r /di; 



(8) 



If the grain size distribution is predominantly modified by 
the three processes considered in this paper, we expect that 
/tot = 1. The deviation of / to t from 1 is an indicator of 
goodness of our assumption that the grain size distribution is 
modified by the three processes. 
2.3 Best fitting parameters 

We search for a set of parameters, (/ grow , /disr), which 
minimizes the square of the difference: 



1 N 
= Jj X! [l°g n synt(a l ) - logn M RN(ai)] 2 



(9) 



where a. b is the grain size sampled by logarithmic bins (i.e., 
log Oi+i — log a, is the same for any i), and N is the number 
of the sampled grain radii (N = 512 in our model, but the 
results are insensitive to N). 

The individual components of processed grain size distri- 
butions [n grow (a, igrow), ^shock(a), and n disr (a, £ d isr)] as 
well as n gi s(a) for silicate and graphite are shown in Fig. QJ 
The MRN size distribution, which should be fitted, is also 
presented. Our fitting procedure is first applied separately 
for silicate and graphite, although we discuss a possibility 
that both species have common values for (/ gIOW , /disr) later 
in Section |4] 

3. Results 

In Table QJ we show the best-fitting values of / grow and 
/disr- We examine t gIOW = 30 and 10 Myr, and tdisi = 3 
and 10 Myr as mentioned in Section |2~T1 We observe that 
/grow = 0.16-0.55 and / diS r = 0.06-0.57 fit the MRN grain 
size distribution. The sum of all the fractions, / tot (Eq. [8]l is 
unity with the maximum difference of 15% (see the column 
of /tot in Table [Tji- 

In order to show how the synthetic grain size distribu- 
tions reproduce the MRN size distribution, we present Fig. [2] 
where we only show Models A and D for the smallest and 
the largest residuals 5 2 . We observe that the best-fitting re- 
sults are fairly consistent with the MRN size distribution. In 
particular, the enhanced and depleted abundances of small 
grains at a < 0.001 fim in naisr and n gs , respectively, cancel 
out very well, especially in Model A. In Model D, the syn- 
thesized size distribution slightly fails to fit the MRN around 
a ~ 0.001-0.002 (im because both n gs and naisr (which are 
used for the fitting) show an excess around this grain radius 



range; thus, the excess around these sizes in the synthesized 
grain size distribution inevitably remains in Model D. How- 
ever, the Milky Way extinction curve is reproduced even by 
Model D within a difference of ~ 10% (see Section FT4l >. 

In order to see the details of the fitting, we show the ratio 
between the synthesized grain size distribution and the MRN 
distribution in Fig. [3] There is a general trend of excess 
around a ~ 0.001-0.003 /im, which is due to grain growth 
(see Fig.Q}. The excess is stronger in Models B and D than 
in Models A and C, which is why the fit is worse in Models 
B and D than Models A and C (Table [Tj- Since the bump 
comes from grain growth, the fit tends to suppress / gr0 w in 
the presence of a strong bump. As a result, / tot is smaller 
in Models B and D than Models A and C (TableQ)- We also 
observe in Fig.[3]that the bump appears at different grain radii 
between Models A/C and B/D because of the difference in 
the duration of grain growth. This bump may disappear 
if coagulation is more efficient than assumed here: more 
efficient coagulation may be realized if i grow 3> 30 Myr 
and/or coagulation also occurs in denser regions (Section 
S3. 

Fig. [3] also indicates that the synthetic grain size distri- 
butions tend to be deficient at the largest grain sizes (a > 
0.1 /im). This is because shattering tends to process large 
grains into small sizes (see Fig. [TJ). The deficiency of large 
grains may be overcome if we include the supply of large 
grains by stellar sources of efficient coagulation as discussed 
in Section l4~2l Because of significant grain growth in Mod- 
els A and C, the deficiency of large grains is recovered by 
grain growth at a ~ 0.01-0.03 ^m for silicate. In Models 
A and C of graphite, shattering causes a dip feature around 
a ~ 0.03 ^m as seen in Fig. [TJ which also appears in Fig. 
[3] In the WIM, where shattering is assumed to occur in this 
paper, grains with a > a few x 10 -2 /zm are accelerated 
up to velocities larger than the shattering threshold by tur- 
bulence. Shattering efficiently destroys small grains because 
of their large surface-to-volume ratios. Thus, the shattering 
efficiency is the largest for the smallest grains that attain a ve- 
locity above the shattering threshold. This is the reason why 
the grains around a ~ 0.03 ^m are particularly destroyed by 
shattering. This dip feature would be smoothed out in reality 
since the grain velocity driven by turbulence has a depen- 
dence on grain charge, gas density, magnetic field, etc., all of 
which have a wide range within a galaxy. 

4. Discussion 

4.1 Derived parameters 

The obtained values of the parameters / gl0 w and /di sr re- 
flect the fraction of individual grain processing mechanisms. 
In other words, these two quantities show the relative impor- 
tance of grain growth and disruption. Note that the efficiency 
of shock destruction is automatically constrained by the bal- 
ance with the mass growth by grain growth (Eq.[6}. Table [TJ 
shows that the best-fitting parameters are not very sensitive 
to t gr ow (duration of grain growth) but that they are sensi- 
tive to tdisf For larger £disr> only a smaller /di sr is necessary 
because the grain size distribution is more modified. As ex- 
pected, /disridisr is less sensitive to idisri note that /disr^disr 
is the mean duration of disruption per processed grain. 

Seeing all the models, we find / gr0 w ~ 0.2-0.6 and 
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Fig. 1. Individual components for synthesized grain size distributions. The thin solid, dashed, and dot-dashed lines represent individual components 
processed by disruption (shattering) for 3 Myr in Panels (a) and (b), and for 10 Myr in Panels (c) and (d), growth for 10 Myr in Panels (a) and (b) and 
for 30 Myr in Panels (c) and (d), and shock, respectively. The thick solid line shows n gj3 (a) = n grow (a, tgrow) + 0.8n s i loc j t (a). The dotted line 
shows the MRN size distribution adopted in this paper. Panels (a) and (c) present silicate while Panels (b) and (d) show graphite. 



/disr ~ 0.06-0.6 (or /disiidisr ~ 0.6-1.7 Myr). As men- 
tioned in Section |22] if silicate and graphite are well mixed 
in the ISM, they are expected to have the common values for 
/grow, /shock, and /disr- In this sense, Models C and D work 
better than Models A and B. In summary, 20-60% of pro- 
cessing occurs in dense clouds (i.e., grain growth), while a 
processed dust grain experiences disruption for ~ 1 Myr on 
average (or disruption accounts for 6-60% of processing). 
From the equilibrium constraint of the total dust mass (Eq. 
|6j, the fraction of shock destruction to all the processing is 



/shock — 0.8/g 



0.1-0.4. 



The sum of all the fractions, / to t (Eq. ©, is unity with 
the maximum deviation of 15%. In other words, we cannot 
reject other processing mechanisms, which could contribute 
to the grain processing with < 15%. 

We have shown that the grain size distributions after (i) 
grain growth, (ii) shock destruction, and (iii) grain disruption 
can synthesize the MRN size distribution. It is also likely 
that we can say the opposite; that is, to realize the MRN size 
distribution, those three processes are crucial. Without (i) 
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Table 1. Models. 



Name 


species 


£grow 

(Myr) 


f'disr 

(Myr) 


/grow 


/disr 


5 2 
(lO" 3 ) 


/tot 


A 


silicate 


30 


3 


0.26 


0.46 


3.7 


0.93 




graphite 






0.43 


0.24 


4.2 


1.01 


B 


silicate 


10 


3 


0.16 


0.57 


16 


0.86 




graphite 






0.40 


0.20 


13 


0.92 


C 


silicate 


30 


10 


0.42 


0.16 


6.4 


0.92 




graphite 






0.55 


0.073 


7.3 


1.06 


D 


silicate 


10 


10 


0.40 


0.15 


33 


0.88 




graphite 






0.50 


0.057 


15 


0.96 



Note: /shock = 0.8/g r ow from Eq. (|6j 
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Fig. 2. Best-fitting synthetic grain size distributions to the MRN size distribution. Panels (a) and (b) show silicate and graphite, respectively. We only 
show two models (A and D; solid and dot-dashed lines, respectively) for the smallest and largest residuals (<5 2 ) among the four models. The dotted line 
shows the MRN size distribution. 



the grain mass just decreases; without (ii) the grain mass just 
increases; without (iii) there is no mechanism that produce 
the large abundance of small grains. 
4.2 Stellar sources? 

In this paper, we did not consider dust supply from stars, 
because the dust mass in the Milky Way is governed by 
the equilibrium between grain growth in molecular clouds 
and grain destruction by SN shocks (e.g. Inoue, 2011). 
Dust grains supplied from stars may be biased to large 
sizes. Production of large grains from asymptotic giant 
branch (AGB) stars is indicated observationally (Groenewe- 
gen, 1997; Gauger et al, 1998; Hofner, 2008; Mattsson & 
Hofner, 2011). The dust ejected from SNe is also biased 
to large grain sizes because small grains are selectively de- 
stroyed by the shocked region within the SNe (Bianchi and 
Schneider, 2007; Nozawa et al, 2007). Coagulation asso- 
ciated with star formation is also a source of large grains if 



coagulated grains in circumstellar environments are some- 
how ejected into the ISM. For this possibility, Hirashita and 
Omukai (2009) have shown that dust grains can grow up to 
micron sizes by coagulation in star formation (see also Ormel 
et al, 2009). As mentioned in Section [3] efficient coagu- 
lation may also solve the bump problem around a ~ 0.001- 
0.003 /zm. These possible sources of large grains may be 
worth including in dust evolution models in the future. 
4.3 Fitting under other constraints 

In Section 12.21 we adopted the balance between the dust 
mass growth by accretion and the dust mass loss by shock 
destruction (Eq.|6]l as a constraint. Although this constraint is 
reasonable for the dust content in the Milky Way (e.g., Inoue, 
2011), it may be useful to apply other constraints without 
using Eq. ©, to see how the best-fitting parameters have 
been controlled by Eq. ©. 

First we try to fit the MRN size distribution with the three 
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Fig. 3. The ratio of the synthetic grain size distribution to the MRN size distribution. Panels (a) and (b) show silicate and graphite, respectively. The solid, 
dotted, dashed, and dot-dashed lines represent Models A, B, C, and D, respectively. 



Table 2. Models with /tot = 1. 



Name 


species 


£grow 

(Myr) 


^disr 

(Myr) 


/grow 


/disr 


S 2 

(un 3 ) 


/tot 


A 


silicate 


30 


3 


0.19 


0.45 


2.1 






graphite 






0.32 


0.24 


3.5 




B 


silicate 


10 


3 


0.065 


0.43 


6.8 






graphite 






0.19 


0.22 


6.2 




C 


silicate 


30 


10 


0.34 


0.16 


4.5 






graphite 






0.53 


0.074 


7.8 




D 


silicate 


10 


10 


0.22 


0.14 


20 






graphite 






0.33 


0.060 


10 





Note: /shock 



/shock- 



Table 3. Models with three free parameters. 



Name 


species 


£grow 

(Myr) 


Edisr 

(Myr) 


/grow 


/shock 


/disr 


8 2 
(10- 3 ) 


/tot 


A 


silicate 


30 


3 


0.13 


0.62 


0.45 


0.98 


1.2 




graphite 






0.20 


0.76 


0.25 


1.8 


1.2 


B 


silicate 


10 


3 


0.038 


0.83 


0.47 


2.4 


1.3 




graphite 






0.11 


0.92 


0.23 


2.1 


1.3 


C 


silicate 


30 


10 


0.27 


0.80 


0.16 


3.0 


1.2 




graphite 






0.42 


0.69 


0.074 


6.4 


1.2 


D 


silicate 


10 


10 


0.16 


1 


0.15 


10 


1.3 




graphite 






0.24 


0.99 


0.061 


5.8 


1.3 



components under the condition that the sum of all the frac- 
tions is unity: / grow + / shoc k + /disr = 1. The results with 
this fitting are shown in Table [2] The best-fitting values of 
/disr vary from those in Table [TJ within a difference of 10% 
except for Model B of silicate (25% less). However, the best- 
fitting values of / gr0 w is broadly 1/2-2/3 of those in Table [TJ 
As a result, Eq. (0 is not satisfied, and the total dust mass 
decreases. 

Next, we perform fitting to the MRN size distribution with 



the parameters f gTOW , /shock, and /disr free- The results are 
shown in Table [3] Again, the values of /di sr differ by only 
< 10% except for Model B of silicate (18% less). However, 
/grow is only ~ 1/3-2/3 of the values in Table [TJ and /shock 
is made large to compensate for the decreased f gIow . This 
means that the fitting is practically dominated by the balance 
between the decreased small grains in shock destruction and 
the increased small grains in disruption (shattering). Because 
of the dominance of /shock, Eq. (O is not satisfied, and the 
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Fig. 4. Upper panel: Extinction curves (extinction per hydrogen nucleus as 
a function of wavelength) calculated for the components used for the fit- 
ting to the grain size distribution. These components are shown in Fig.fT] 
The thin solid, dashed, and dot-dashed lines represent individual compo- 
nents processed for the following processes: disruption (shattering) for 3 
Myr in Panel (a) and for 10 Myr in Panel (b), growth for 10 Myr in Panel 
(a) and for 30 Myr in Panel (b), and shock, respectively. The thick solid 
line shows the extinction curve for n g , s (a)/1.8 (divided by 1.8 because 
the component "g,s" contains the grain mass 1.8 times as much as the 
MRN). The dotted line presents the extinction curve for the MRN size 
distribution. The points show the observed Milky Way extinction curve 
taken from Pei (1992). Lower panel: Ratio of extinction curves to the ex- 
tinction curve of the MRN size distribution. The line species in the lower 
panel correspond to those in the upper panel. 



total dust mass decreases. 
4.4 Extinction curve 

The MRN grain size distribution is originally derived from 
the Milky Way extinction curve. Therefore, in order to check 
if our fitting by synthetic grain size distributions is successful 
or not, it is useful to calculate extinction curves. 

Extinction curves are calculated by using the same optical 
properties of silicate and graphite as those in Hirashita and 
Yan (2009). The grain extinction cross section as a function 
of wavelength and grain size is derived from the Mie theory, 
and is weighted for the grain size distribution per hydrogen 
nucleus to obtain the extinction curve per unit hydrogen nu- 
cleus (denoted as A\/Nh). The abundances of silicate and 



graphite relative to hydrogen nuclei are already inherent in 
the models through the abundances of Si and C and £ (Sec- 
tion [2j. 

First we show the extinction curves of the individual com- 
ponents, which are used to fit the MRN size distribution, in 
Fig- SI Grain growth does not make the extinction curve flat- 
ter in spite of the increase of the mean grain size. The rea- 
son is already explained in Hirashita (2012): Accretion pre- 
dominantly occurs at the smallest sizes. Since the extinction 
at short wavelengths is more sensitive to the increase of the 
mass of small grains than that at long wavelengths, the ex- 
tinction curve becomes rather steeper. Although coagulation 
flattens the extinction curve, the flattening due to coagulation 
does not overwhelm the steepening due to the above effect of 
accretion. 

Shock destruction makes the extinction curve flatter be- 
cause small grains are more easily destroyed than large 
grains. Grain disruption steepens the extinction curve be- 
cause of the production of a large number of small grains. 
The 0.22 /im bump created by small graphite grains in this 
model becomes also prominent by grain disruption. We 
also show the extinction curve for the grain size distribution 
%,s(a) (The component n gs has a total dust mass 1.8 times 
as large as the initial value. To see the difference of the ex- 
tinction curve, it would be helpful to compare under the same 
dust mass, so n gs /1.8 is compared with the MRN in Fig. [4]) 

In Fig. |5j we show the extinction curves calculated for 
Models A-D. First of all, we confirm that the MRN size 
distribution reproduces the observed extinction curve (some 
small deviations can be fitted further if we adopt a more de- 
tailed functional form of the grain size distribution, which is 
beyond the scope of this paper; see Weingartner and Draine, 
2001, for a detailed fitting). Comparing the extinction curve 
for the MRN size distribution and those for Models A-D, we 
observe that the extinction curve is reproduced within a dif- 
ference of ~ 10%. Models A and C are successful, while 
Models B and D systematically underproduce the MRN ex- 
tinction curve (although the difference is small). The under- 
prediction by ~ 10% in Models B and D occurs because /tot 
is - 0.9. 

In the above, we adopted different values for / gr0 w and 
/disr between silicate and graphite. As mentioned in Sec- 
tion 12.21 if both species are well mixed in the ISM, they 
would have common values of these parameters. In Fig. [6] 
we show the extinction curve by taking the average of the 
values for silicate and graphite (for example, / gr0 w = 0.35 
and /disr = 0.35 for Model A). We find that the differ- 
ence between Figs. [5] and [6] is small. Therefore, the mean 
values work to reproduce the Milky Way extinction curves. 



The mean values are in the range of / g 



0.3-0.5 and 



/disr = 0.1-0.4. We conclude that the synthetic grain 
size distributions with these parameter ranges reproduce the 
Milky Way extinction curve. 

5. Conclusion 

In our previous papers (Nozawa et al., 2006; Hirashita 
and Yan, 2009; Hirashita, 2012), we showed that dust grains 
are quickly processed by shock destruction, disruption, and 
grain growth. In this paper, thus, we have examined if the 
MRN grain size distribution, which is believed to represent 
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Fig. 5. Extinction curves (extinction per hydrogen nucleus as a function 
of wavelength) calculated for Models A-D (upper panel). The ratio to 
the extinction curve for the MRN size distribution is also shown (lower 
panel). The solid, thick dotted, dashed, and dot-dashed lines represent 
Models A, B, C, and D, respectively. The thin dotted line shows the 
extinction curve for the MRN size distribution. The points show the 
observed Milky Way extinction curve taken from Pei (1992). 
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Fig. 6. Same as Fig. [5] but we use the mean values between silicate and 
graphite for / grow and / disr . 



the grain size distribution in the Milky Way, can be repro- 
duced by the processed grain size distributions. We "synthe- 
sized" the grain size distribution by summing the processed 
grain size distributions under the condition that the decrease 
of dust mass by shock destruction is compensated by grain 
growth. We have found that the synthetic grain size dis- 
tribution can reproduce the MRN grain size distribution in 
the sense that the deficiency of small grains by grain growth 
and shock destruction can be compensated by the produc- 
tion of small grains by disruption. The values of the fit- 
ting parameters indicate that, among the processed grains, 
30-50% is growing in dense medium, 20-40% is being de- 
stroyed by shocks in diffuse medium, and 10-40% is being 
shattered in diffuse medium (the percentage shows the rel- 



ative importance of each process). The extinction curves 
calculated by the synthesized grain size distributions repro- 
duce the observed Milky Way extinction curve within a dif- 
ference of ~ 10%. This means that our idea of synthesizing 
the grain size distribution based on major processing mecha- 
nisms (i.e., grain growth, shock destruction, and disruption) 
is promising as a general method to "reconstruct" the extinc- 
tion curve. 
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